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A nonlinear model of spin wave excitation using a point contact in a thin ferromagnetic film 
is introduced. Large-amplitude magnetic solitary waves are computed, which help explain recent 
spin-torque experiments. Numerical simulations of the fully nonlinear model predict excitation 
frequencies in excess of 0.2 THz for contact diameters smaller than 6 nm. Simulations also predict 
a saturation and red shift of the frequency at currents large enough to invert the magnetization 
under the point contact. The theory is approximated by a cubic complex Ginzburg-Landau type 
equation. The mode's nonlinear frequency shift is found by use of perturbation techniques, whose 
results agree with those of direct numerical simulations. 
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Generation of spin waves at microwave frequencies due 
to spin-momentum transfer (SMT) in a thin magnetic 
multilayer is paving the way to a new frontier in magne- 
todynamics. Initial predictions, in the small amplitude 
limit, made by Slonczewski 0] were recently validated 
experimentally El In refs 0, 01 , microwaves were 
generated by the application of a dc current through a 
nanocontact into a trilayer of alternating ferromagnetic 
and nonmagnetic materials (see fig Several mod- 
els H, IS 01 have been suggested to explain this phe- 
nomenon. More recent work has neglected spatial vari- 
ation 0, 0, lld| . In this study, a theory of nonlinear 
spatially nonuniform magnetodynamics is developed and 
analyzed. Microwave emission from point contacts has 
potential applications in communications and fast, high- 
density magnetic data storage. 




FIG. 1: Magnetic multilayer with magnetization M precessing 
in the free layer with frequency il). The magnetization of the 
lower ferromagnetic layer is fixed with orientation rfiF- A 
nonmagnetic spacer is sandwiched between two ferromagnetic 
layers. Typical parameters are: r* = 20 nm and S — 5 nm. 

The system depicted in fig. ^ nas been shown exper- 
imentally to give rise to high-frequency microwaves in- 
ferred to be caused by steady-state precession of magne- 
tization in the ferromagnetic free layer (NigoFe2o) [m- An 
external magnetic field Sq is oriented perpendicular to 
the film plane. When a dc current, I, is applied through a 



point contact into the multilayer CogoFeio/Cu/NigoFe2o, 
the conduction electrons' spins induce a torque on the 
magnetization in the free layer due to conservation of 
angular momentum P, 0] . This SMT torque opposes the 
intrinsic damping of the free layer and produces complex 
dynamics. Slonczewski's linear theory of spin wave exci- 
tation for point contacts on magnetic multilayers predicts 
a threshold current for the excitation of steady precession 
of magnetization However, it does not accurately 

predict the high frequencies seen in large amplitude ex- 
periments nor the dependence of the frequency on cur- 
rent. 

Understanding the SMT effect in point contacts with 
a perpendicular geometry is of great interest as it is a 
microwave source that yields the narrowest line widths 
and highest output powers 0, 0] when compared to 
other SMT geometries, such as nanopillars Q and in- 
plane fields |13| . In addition, this configuration is most 
amenable to analytical study given the high symmetry of 
the problem. Insights derived from this simpler geome- 
try may then be employed in solving the more advanced 
problems of in-plane or oblique angled fields. 

In this Letter, beginning with a nonlinear vectorial 
equation, we derive a new complex cubic Ginzburg- 
Landau (CGL) equation for weakly nonlinear excita- 
tions. The dipole coupling and exchange nonlincaritics 
are shown to be responsible for the higher frequencies. 
Steady-state precessing modes for the nonlinear vecto- 
rial model are found via numerical integration. The cal- 
culated frequencies lie in the range of experimental ob- 
servations jj] . Using this model, much higher frequencies 
(0.2 THz) are predicted for systems with contact radii on 
the order of r* = 3 nm. A surprising prediction is the 
saturation and red shift of frequency for currents large 
enough to invert the magnetization under the point con- 
tact. This is a direct consequence of the exchange term 
included in our model. 

The magnetization M = (M x , M y , M z ) in the free layer 



2 



can be described by the vectorial equation 

^ = - | 7 | /i M x R^--^—M x (M x tf eff ) ■ 
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Landau-Lifshitz damping 
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Slonczcwski SMT torque 

where the precessional and damping terms are driven 
by the effective magnetic field i? e ff = Hqz — M z z + 
\i\JaM n ^ 2 This field consists of the applied magnetic 
field, Hqz, the demagnetizing field due to axial dipole 
coupling, —M z z, and the exchange field proportional to 
V 2 M 14] . The Landau-Lifshitz form of damping is a 
commonly used phenomenological term that drives the 
magnetization to align with the total field. The Slon- 
czewski SMT torque term is derived in and assumes a 
large, bulk ferromagnetic layer with fixed magnetization 
direction i. Relevant parameters are the gyromagnetic 
ratio (7 = g/is/fi.), the spectroscopic splitting factor (g), 
the Bohr magneton (/is), Planck's constant divided by 
2ir (h), the free space permeability (po), the exchange 
coupling parameter (-D), the magnetization volume den- 
sity at saturation (M s ), and the transverse relaxation 
time (T2); (3(x) is a driving term depending on the cur- 
rent discussed below. 

We assume zero temperature and negligible crystalline 
anisotropy in order to gain insight into the basic mode 
structure facilitated by the fundamental physics associ- 
ated with SMT in point contacts. Such effects can be 
accounted for by modification of eq JQ) to include a ran- 
dom fluctuation term and an effective anisotropy field. 
Also, since the thickness of the free layer is small, S m 5 
nm, the in-plane components of the dipole field contri- 
bution to H c g can be neglected |l5j . 

The magnetization is assumed to have rotational sym- 
metry: the spatial variation of the magnetization de- 
pends solely on the distance r from the center of the 
point contact. Then, the current term, f3, is given by 
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where e is the SMT efficiency, 6 is the thickness of the 
free layer, e is the charge of an electron, r — {x 2 +y 2 ) 1 / 2 , 
and $(r* — r) is the Heaviside step function defining the 
point contact to be a circle of radius r* . The efficiency 
e is in principle a complicated function of many param- 
eters, including (but not limited to) microscopic details 
of the interfaces in question 0, [n| . Many of the approx- 
imations used in the calculation of e ignore the lateral 
geometry of the actual experiments, choosing instead to 
use a more tractable ID calculation. Since e is difficult 
to approximate, let alone calculate, and is not the sub- 
ject of this paper, we will treat e as a constant fitting 
parameter. 



By taking the dot product of equation with M, one 
sees that the magnetization of the free layer is locally 
conserved for all times, i.e., \M\ — M s — constant, or 

M 2 S = M*(r,T)+MZ(r,T)+M*(r,T), r > 0, r > 0. (2) 



Consider the standard normalization, m = M/M s — 
(m x ,m y ,m z ), t = u! M r, p = r/ l CXl where aj M = TMoMs 
and the exchange length is ^ cx = y/ D /jp M s h. For the 
analysis, it is convenient to encode the transverse com- 
ponents of m, m± = (m x ,m y ), in the complex quantity 
m = m x + irriy. The axial component of the magneti- 
zation is treated as a perturbation from the equilibrium 
solution m = z, so m z = 1 — /, where / is the normalized 
perturbation. Evaluating the vector cross products in 
leads to 



= toV 2 / + (1 - / + ia)V 2 m + (1 - / - h)m 
+i{j$(p* -p)- a(h - (1 - /))}(1 - f)m 
-ia\ |V/| 2 + |Vr7j| 2 \m 



^ = Im(m*\7 2 m) + {j<5>(p*-p)-a(h-(l-f))} \mf 



-a{v 2 /-(|Vm| 2 + |V/| 2 )(l-/)}, 



(3) 



where V 2 = + -■§-• The dimensionless parameters 
for applied field, damping, and current are respectively 
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The coupled system © is equivalent to the model QJ. 
Since the patterned multilayer mesa used in experiments 
is two orders of magnitude wider than the point contact 
fl, we assume that the multilayer has infinite extent in 
the xy plane, hence < p < 00. 

The nominal parameter values used in this paper are 
r* = 20 nm, e = 0.26, D = 4mcV-nm 2 , M s = 640 kA/m, 
6 = 5 nm, 7=1.85T0 n Hz/T, w M = 23.68 GHz, Z cx = 6.40 
nm, and a = 0.0112. We show that by using only e as a 
fitting parameter, the theoretical data match experiment 
well. Except where noted, the parameter values h=l.l 
and =3.12 are used. 

In order to develop a perturbation scheme, we intro- 
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duce a small parameter, ^-<Cl, and a rescaling of equa- 
tion |J3J|, m=arh. The parameter a represents the magni- 
tude scale of the transverse magnetization at the center of 
the point contact. Using the approximation, / w ^-|m| 2 
for the saturation condition, substituting this approxi- 
mation into eq ©, and keeping only first order terms in 
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a 2 , we find that m satisfies a CGL-type equation, 



dm 
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dipolc coupling 



j& — a(h — 2) |m| 2 rh, 



mV z m* + 2(1 + ia) |Vm| 



(4) 



nonlinear exchange 



Seeking a steady-state solution of the form m(p, t) — 
e ItAjt 4>(p), Cj = h—l—u) yields a nonlinear eigenvalue prob- 
lem with boundary conditions 0(0) = 1, dtfi/dp(0) = 0, 
linip^oo <ft(p) = 0. The corresponding linear eigenvalue 
problem was solved by Slonczewski with two real 
eigenvalues u> and j. This suggests the asymptotic ex- 
pansions 



(5) 



LJ = W (0) + — LJ 



,(1) 



i=i (0) + yj 
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Slonczewski's linear res ult is summarized below. 
Linearized solution [11]: The leading order mode 
'(p) in the expansion 0) is 



Jo(hp) < p < p* 
p) p*< p 



, = ( - L ,m +l ( j (°)- a (h-l)) )1 /2 k _ ( -^-i a (h-l) )1/2 

where the linear eigenvalues ui^ and are determined 
by the real and imaginary parts of 



faHQ 2 \k p*)Ji(kip*) = koHf (k o p*)J (kip*). (6) 

Equation © is solved numerically using a nonlinear 
root finder. A doubly-discrete set of eigenvalues for fre- 
quency and current is found. We consider only the lowest 
eigenvalues o/°J andj^°\ where is the critical current 
for linear spin wave excitations 11]. 

The predicted frequency of the ground state mode, 
Cj ~ h — 1 — uj(°\ is valid only for very small amplitudes 
a, whereas experiments produce much higher frequencies 
|4| . The nonlinearity in (@J gives rise to a frequency shift 
which, through the use of the Poincare method 0], helps 
explain this disparity. 

Nonlinear solution: The first order nonlinear fre- 
quency and critical current shifts a/ 1 -* and are de- 
termined by the linear system 



. rr(2), 
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where Ik depend only on the leading- order solution: 
h = h-^((4 2 \k p r )) 2 + (H[ 2 \k p r )) 2 



-|^| 2 + ^V 2 ^* + 2|V0(°)| 2 j0(°) pdp 

a(2 - h)\^\ 2 - 2a|V0W| 2 |^ 2 (O dp 
+ j (0) T* | (o)|2 (o) 2 pcip _ 



The above perturbation method provides a means for 
calculating the nonlinear frequency and current shifts, 
a/ 1 ) and j", of a mode solution to 10} and helps explain 
experimental results 0. The frequency and critical cur- 
rent have a parabolic dependence on the transverse am- 
plitude a 

fi-fc-l-O^-^W, jr.jW+^jW. ( 7 ) 

We note that using equations J7J), the relation between Cj 
and j is linear because both scale as a 2 . 

The first-order nonlinear contributions from equation 
Q due to exchange and dipole coupling were compared. 
Using equations Q), the frequency and current depen- 
dence as parameterized by a were calculated when the ex- 
change and dipole coupling nonlinearitics were removed 
from the calculation. We find that if nonlinear exchange 
terms are ignored, the relative error in the slope dd/dj is 
approximately 10% for r* = 20 nm, whereas for < 12 
nm the relative error is greater than 25 %. 

Experiments show that for large currents, and hence 
larger amplitudes, the frequency/current relationship sig- 
nificantly deviates from linearity. In order to account for 
this, equation (|2J) was solved numerically using the pro- 
jection method where one renormalizes the magnetiza- 
tion after each time step in order to preserve the local 
conservation of M @ (see [13 )■ The steady states of 
this parametrically forced/damped system are attractors. 
The initial condition m(p, 0) = a4>^\p) was evolved in 
time over a large domain (0 < p < 100) until the solu- 
tion relaxed to a steady state. The mode frequency was 
numerically recovered from the phase of the precessing 
transverse component as Cj = ^ arg[m(0, t)}, t — > oc. 

Typical modes, their frequencies, and currents are 
shown in fig. [21 The transverse mode envelope ± \m\ = 
±(m 2 + m 2 ) 1 / 2 and the y-component of magnetization 
m v show slow decay, which means that spin waves are 
continually radiating away from the point contact, and 
hence energy from the contact region is lost into the sur- 
rounding medium. 

The plot of frequency versus current in fig. |2| shows a 
comparison between numerical simulation of J2J, pertur- 
bation theory J7J) , and experiment Q . For current excita- 
tions larger than j max , the solitary- wave mode structure 
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FIG. 2: Spatial dependence of fully nonlinear modes with 
excitation current and frequency. Top: transverse envelopes 
±|m| = ±^/ml + rriy and associated spin waves m y . Bottom: 
axial components m z . For j > j ma x, the transverse mode 
magnitude is no longer monotonically decreasing since the 
axial component of magnetization is negative (m z < 0) in the 
point contact. Frequency and current values for each mode 
are shown as circles in fig. [3 



becomes highly nonlinear (compare with fig. [2J. There 
is a maximum frequency and corresponding current for a 
given physical system (jmax,&max)- As the dashed mode 
in fig. 0shows, frequency saturation corresponds to mag- 
netization precession along the equator near the center of 
the point contact (m z « 0). For currents less than j max , 
the transverse mode magnitude has a monotonically de- 
creasing dependence on r. For larger currents, the mode 
is not monotonic and develops interesting structure in- 
side the point contact. 

Fig. |2l shows that, as the point contact size is in- 
creased, the slope of frequency versus current decreases 
because the slope depends on SMT torque (inversely pro- 
portional to rl) and damping due to spin wave genera- 
tion (inversely proportional to 7'* in the linear approxi- 
mation ^3)- The comparison of the simulated results (*) 
with experiment is remarkable when the SMT efficiency 
is taken to be e = 0.8. As the experimental values of 
the frequency level out, full saturation is being reached. 
This is the first model to qualitatively capture this type 
of behavior. 

Nominal values for e are on the order of 0.25 for most 
ferromagnetic metals ^8|. The fact that e > 0.25 is used 
to fit the experimental data suggests that further refine- 
ment of the model is desirable. 

The dependence of frequency on applied magnetic field 
H is depicted in fig. g| With the fitted SMT efficiency 
parameter, the theory compares well with experiment. 
The mode structure is highly nonlinear because the trans- 
verse amplitude is |m(0,f)| ~ 0.98, which precludes the 
use of eq (0). However, a full numerical solution is re- 
quired for only one particular value of h. Since duj / dh ss 
constant for a«l. once the parameters a, j, and p* 
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FIG. 3: Precession frequency as a function of current with 
comparison between fully nonlinear modes (direct numerical 
simulation of eq the perturbation result (J, and exper- 
iments conducted using r, — 20 nm point contacts Q. Large 
circles represent the frequency and current of modes plotted 
in fig. [5] The straight lines were generated by perturbation 
theory (eq (J7J) for the contact radii denoted. Note the excel- 
lent agreement between numerical simulation, perturbation 
theory, and experiment for small currents. The parameters 
have nominal values except e = 0.8. 
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FIG. 4: Frequency as a function of applied field h. Compar- 
ison between fully nonlinear modes (direct numerical simula- 
tion of eq JjU) and experiments 4]. Parameters are I = 10 
mA and fitted e = 0.5. Inset: maximum frequency (*) and 
current (o) versus contact radius. 



are chosen, the y-intercept for the approximate linear re- 
lationship between lj and h is fixed. The slope of the 
line does not change. The theory predicts the slope very 
accurately. In fig. ^ inset, much higher frequencies are 
predicted by shrinking the size of the point contact. Us- 
ing the method stated above, we predict that frequencies 
in the 0.2 THz range are attainable for a point contact 
with a radius of 3 nm. We find that as the point contact 
is made smaller, the current decreases linearly for the 
same reasons stated earlier in regard to the Cj(j) slope, 
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while the frequency increases in proportion to 1/rJ; due 
to exchange- mode spin- wave dispersion. 

In summary, a new model of magnetic excitations in 
point contact structures that includes nonlinearities due 
to both exchange and dipole coupling is introduced. Spa- 
tially dependent steady-state modes are calculated. This 
represents a large step forward in the explanation of the 
phenomena observed in recent experiments. 
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